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When a granular mixture involving grains of different sizes is shaken, sheared, mixed, or left to 
flow, grains tend to separate by sizes in a process known as size segregation. In this study, we 
explore the size segregation mechanism in granular chute flows in terms of pressure distribution and 
granular micro-structure. Therefore, 2D discrete numerical simulations of bi-disperse granular chute 
flows are systematically analysed. Based on the theoretical models by Gray & Thornton [33] and 
Hill &; Tan [31], we explore the stress partition in the phases of small and large grains, discriminating 
between contact stresses and kinetic stresses. Our results support both gravity-induced and shear- 
gradient-induced segregation mechanisms. However, we show that the contact stress partition is 
extremely sensitive to the definition of the partial stress tensors, and more specifically, to the way 
mixed contacts (i.e. involving a small grain and a large grain) are handled, making conclusions on 
gravity-induced segregation uncertain. By contrast, the computation of the partial kinetic stress 
tensors is robust. Kinetic pressure partition exhibits a deviation from continuum mixture theory of 
a significantly larger amplitude than contact pressure, and display a clear dependence on the flow 
dynamics. Finally, using a simple approximation for the contact partial stress tensors, we investigate 
how contact stress partition relates with the flow micro-structure, and suggest that the latter may 
provide an interesting proxy for studying gravity-induced segregation. 

PACS numbers: 45.70.-n, 05.65.+b 


I. INTRODUCTION 


When a granular mixture involving grains of different 
sizes is shaken, sheared, mixed, or left to flow, grains 
tend to separate by sizes in a process known as size seg¬ 
regation. In the simple case of a free surface flow under 
gravity, larger grains rise to the surface while the smaller 
grains sink at the bottom of the flow. Segregation oc¬ 
curs in a wide variety of contexts - either industrial or 
natural inn - for which it represents an engineering and 
scientific challenge. Kinetic theory offers a well estab¬ 
lished framework to describe segregation in dilute gran¬ 
ular systems 0 0 ; however its applicability is limited 
in the case of dense flows, where grains interact through 
long-lasting contacts rather than binary collisions. Since 
the seminal work by Savage and Lun [9] and the intro¬ 
duction of the random kinetic sieve , much progress has 
been achieved in understanding segregation mechanisms, 
both experimentally |TQH22l and using numerical simula¬ 
tions [23H32] . In this respect, because they allow for very 
well controlled configurations and give access to the inner 
structure of the flow, discrete numerical simulations have 
proven a fruitful tool, and provide significant insight to 
inform continuum modeling [4J EU GSHSS] • 

One of the challenge posed by granular size segregation 
lies in the task of identifying the actual driving mech¬ 
anisms [27] EH [331 [37H39]. Experimental evidence for 


the lift force that drives an intruder upwards in a quasi¬ 
static shear flow was recently given and shown to scale 
like Archimedes force m ■ It seems reasonable to sup¬ 
pose that such lift forces are at play in size-disperse dense 
granular flows, although the dependences on the flow 
composition and dynamics have to be established. In the 
absence of an explicit generic formulation for the segre¬ 
gation force, a fruitful hypothesis was proposed by Gray 
& Thornton [33]: the smaller grains, while percolating 
through the network of larger grains, are screened from 
the average stress and thus carry less of the lithostatic 
pressure than larger grains do. This hypothetical mech¬ 
anism translates mathematically into a departure from 
classical continuum mixture theory; it allows for the re¬ 
covery of important features of segregation in depth av¬ 
eraged equations [33] . It has been since the subject of 
various improvements [291 EDI IMl [38]. Yet the basic in¬ 
gredient, namely the fact that the partial pressure in the 
phase of the smaller grains is less than simply propor¬ 
tional to their volume occupation, has to be unambigu¬ 
ously established. While Weinhart et al do observe such 
asymmetry in discrete numerical simulations [39] , the re¬ 
sults obtained by Fan & Hill and Hill & Tan are inconsis¬ 
tent with this hypothesis mum]. More recently, it was 
proposed by Fan & Hill that in addition to gravity - or 
in place of it - kinetic stresses play a crucial role in size 
segregation m They identify the gradients of velocity 
fluctuations as the driving mechanism [271 EH- 
Based on the theoretical models by Gray & Thornton 
[33] and Hill & Tan [31], we explore the stress partition 
in the phases of small and large grains, discriminating 
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FIG. 1: Example bi-disperse granular flow with volume frac¬ 
tion of large beads = 0.3 in the initial state (t — 0) 
and after steady-state is reached and segregation occurred 
(t = 150s); the slope is 0 — 23°. The ratio between the diam¬ 
eter of large and small grains is d^/ds — 2 (ds — 0.005 m). 
(color online) 


between contact stress and kinetic stress, in numerical 
bi-disperse granular flows. The numerical method and 
simulation are presented in section [ll| Contact stresses 
and their partition between the phases of small and large 
grains are analysed in section m Section [IV] explores 
the partition of kinetic stresses and its sensitivity to the 
flow dynamics. The relation between contact stresses and 
granular micro-structure is established in section [V| We 


finally discuss the results in section VI 


II. THE NUMERICAL FLOWS 


solid friction: locally, the normal and tangential contact 
forces satisfy f t < /i/ n , where p is the coefficient of 
friction at contact. Moreover, a coefficient of energy 
restitution e sets the amount of energy dissipated by 
collisions. The numerical values of p and e control 
the effective frictional properties of the flow in a given 
configuration. They are strictly the same for all contacts 
between large and small grains. In this work, we are 
not interested in understanding how they may affect the 
segregation process, hence, their value was set to p = 0.5 
and e = 0.25, and was not varied. 

We denote dz, and ds the mean diameter of large 
and small grains respectively. To prevent geometri¬ 
cal ordering, likely to happen for strictly mono-sized 
packings, both large and small grains have diameters 
uniformly distributed around their mean value so that 

/ imax_ imin \ / imax_ imin \ 

- - dL - - = - s ds s - = 0.08 (a discussion on the 
influence of this value on the segregation process can 
be found in 22). The ratio d^jds was not varied: 
dhl^s = 2. 

Periodic boundary conditions were implemented to en¬ 
sure long flow durations; the width of the simulation cell 
is 35 dL- The basal boundary is made of a row of fixed 
beads of diameter dz,. In the initial state, a layer of large 
beads is overlaid by a layer of small beads (Figure [TJ. 
This is achieved by random deposition under gravity. 
We denote <f> L the volume fraction of large beads at the 
flow scale, i.e. the ratio of the volume of large beads to 
the total volume of grains: <f> L = Vl/(Vs + Vl). The 
volume fraction of small beads is = (1 — <f> L ). In 
the simulations, <f> L was varied between 0.06 and 0.90. 
In addition, we define the local volume fraction of large 
and small beads (f) L and (j ) S . The height of the granular 
bed in the initial state is Hq; irrespective of <f> L , Hq was 
kept constant and equal to H ~ 60ds- The slope of the 
granular bed 0 was varied between 21° and 26°, allowing 
different flow velocities. The numerical values used for 
the simulations are the following: ds = 5.10 -3 m, p = 1 
kg.m -2 , g — 9.8 m.s -2 . 


III. CONTACT STRESSES AND 
GRAVITY-DRIVEN SEGREGATION 


Contact stresses are transmitted at contacts between 
grains during either short collisions or enduring contacts; 
its micro-mechanical expression is the following 02t 


The numerical systems consist of two-dimensional 
granular beds of grains of two sizes allowed to flow under 
gravity on a fixed plane of slope 0 made of grains of the 
larger size. As the flow develops, segregation occurs, 
and the larger grains - placed initially at the bottom 
of the granular layer - rise in the flow, as illustrated in 
Figure [l] This typical chute flow experiment is in its 
principle similar to those reported in [25j |29l |30j |32] . 
This numerical experiment was performed using the 
contact dynamics algorithm mum, assuming perfectly 
rigid grains. The grains interact at contacts through 


<T = ^7 E (!) 

aeN c 

where f a is the force transmitted at the contact <a, £ a is 
the vector joining the centres of mass of the two grains 
involved, N c is the number of contacts over which the 
summation is made, V is the volume over which the stress 
is computed, and 0 is the dyadic product. The eigen val¬ 
ues Ai and A 2 of this stress tensor gives the pressure P: 
P = (Ai + A 2 V 2 . A typical pressure profile is shown in 
Figure ^a for a granular flow with <f> L = 0.30 flowing at 
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an angle of 23° in the steady (i.e. segregated) state, as 
shown in Figure [l] It simply obeys a lithostatic profile 
P(z) = pg cos 0 (Hq — z), where z is the depth (counted 
from the bottom). 

In a bi-disperse granular flows, contacts may involve only 
large grains, only small grains, or one large grain and 
one small grains. In these last mixed cases, a meaningful 
partition of the stress for the computation of the partial 
stress tensors a s and cr L is necessary. Following HQ [39], 
the contribution of the mixed cases is distributed accord¬ 
ing to the size of the grains, namely weighted by a pref¬ 
actor dL/(dL +ds) or ds/(dL + ds) for the phase of large 
and small grains respectively (we shall see that using a 
different partition changes the results dramatically). We 
thus define the partial stress tensors cr L and cr 5 , and the 
corresponding partial pressures P L and P s . In the fol¬ 
lowing however, to allow precise comparison with earlier 
works, we will mostly use the normal stress components 


(jy y and Gyy 


In the framework of classical continuum mixture theory, 
partial stresses are expected to be proportional to the 
mean stress a yy and to the volume fraction of the con¬ 


sidered grain species in the mixture, namely locally in the 
flow: a yy = (\> L a yy and a yy = <j) S a yy . A departure from 
mixture theory as proposed by 33 implies however dif¬ 
ferent stress partition coefficients vp L ^ <p L and ijj s ^ <fi s 
[271 [39] such that ip L + ip s = 1 and 


a ; 


= t°yy, 

= P°yy 


Considering bi-disperse flows with a volume fraction of 
large grains <I> L = 0.06, 0.15, 0.30, 0.45, 0.60 and 0.75, 
and flowing at an angle 0 = 23°, we evaluate the local 
contact stresses a yy and a yy using equation (Jlj) , com¬ 
puted over horizontal layers of width 4 ds, and over time 
intervalles of 1.25 seconds. For the same volumes and 
time intervals, the local volume fraction of large and 
small grains <p L and <fi s are also evaluated. We do no try 
to separate the early stages of the segregation from the 
later stages as we found that it had no visible influence 
on the results shown here after. We can thus plot the 
stress ratio a yy / a y as a function of volume fraction ra¬ 
tio (j) S /(j) L (Figure®}). In order to filter out the extreme 
cases, the best fit is calculated for 0.01 < (j) S /(j) L < 100. 
We find 


<j^ as 

~ (0.9660 ± 0.003)(2) 

which suggests a small asymmetry of the stress partition 
compared to a classical mixture. 

In the original model proposed by Gray & Thornton [33] , 
the stress partition obeys the following law: 

P = + ( 3 ) 

P = 0 S (1 -B<f> L ). (4) 

Using ([ 2 ]), (( 3 J) and (j4j) gives B ~ (1 — 0.9660)/(0^ + 
0.9660(1 — </T)), namely 0.034 < B < 0.035. This value 



FIG. 2: a Exemple of pressure profile in a segregated flow 
with = 0.3 and 9 — 23°; the dotted line shows P — 
pg(Ho — z); b Ratio of the partial contact stresses cry y /(jyy 
as a function of the local ratio of the volume fractions 4> S /4> L 
for flows with 4> L = 0.06, 0.15, 0.30, 0.45, 0.60 and 0.75, and 
slope 0 — 23°; the line shows a yy / a yy — 0.9660 (0 S /0 L ); c 
Proportion of mixed contacts (between one small and one 
large grain) as a function of the local volume fraction of 
large grains for <f> L = 0.06, 0.15, 0.30, 0.45, 0.60 and 0.75, 
and slope 0 = 23°. 


is larger than the value observed for 3D chute flows by 
Weinhart et al (who find B ~ 0.02) [39] . and supports the 
plausibility of gravity-driven segregation in the system 
(note that Hill & Tan find B 0 for 3D rotating drums 

ED). 


It should be noted that the value of B is very sensitive 
to the details of the stress computation. For instance, 
considering the ratio of the partial pressures P s /P L 
rather than the ratio of the partial normal stresses 
Gyy/Gyy gives 0.045 < H < 0.047, namely a larger value 
of B. Far more dramatic is the influence of the stress dis¬ 
tribution at mixed contacts (between a large grain and 
a small grain) between the two phases. In the analysis 
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above, the distribution of the stress transmitted at mixed 
contacts is weighted by the grain size, namely a prefactor 
of 1/3 for the contribution to the phase of small grains 
and a prefactor of 2/3 for the contribution to the phase 
of large grains (since dj, = 2 ds). Changing the prefac¬ 
tors to 0.2 for the phase of small grains and 0.8 for the 
phase of large grains leads to B ~ 0.33, namely a massive 
increase. On the contrary, changing the prefactors to 0.4 
and 0.6 leads to B ~ —0.1, that is an inverse segrega¬ 
tion process. Finally, changing the prefactors to 0.5 and 
0.5 leads to B ~ —0.32, and brings an entirely different 
picture. The strong influence of the way mixed contacts 
are taken into account in the stress computation can be 
explained by their proportion in the mixture, which can 
reach locally 50% depending on the volume fraction of 
large grains, as is visible in Figure [2^. This makes the 
interpretation of contact stress ratios awkward, partic¬ 
ularly if one wants to compare mixtures with different 
grain sizes. In that case, it will be difficult to tell the 
bias introduced by the ponderation by the grains size in 
the partial stresses computation from the influence on 
the mixed stress partition that derives. 

No influence of the slope angle 0 (namely of the flow ve¬ 
locity) on the value of B was observed. 

IV. KINETIC STRESSES AND 
SHEAR-GRADIENT-INDUCED SEGREGATION 

Kinetic stresses are related to the existence of veloc¬ 
ity fluctuations in granular flows, and can be quantified 
through an analogue of the Reynolds stress tensor: 

cr K = / y] rriiSvi <g> 5v it (5) 

ieN p 

where Svi is the fluctuating velocity of the grain i, rrii its 
mass, N p the total number of grains in the volume V over 
which the stress is computed. A typical kinetic pressure 
profile is shown in Figure^ for a steady state granular 
flow with 4 > l = 0.30 flowing at an angle of 23°. The nor¬ 
malisation by pgds allows for quantitative comparison 
with the contact stress profile shown in Figure [2] Ki¬ 
netic stresses are about four orders of magnitude smaller 
than contact stresses, namely seemingly negligible. How¬ 
ever, as stressed by EEJEH], the difference of behaviour 
between the phase of small grain and the phase of large 
grain is more striking than for contact stresses. Com¬ 
puting the partial kinetic stress tensors cr K,s and cr^ ,L is 
straightforward as it implies a summation on the grains 
and not on the contacts. Plotting the kinetic stress ratio 
a yy S l^yui as a f unc ti° n of volume fraction ratio (j) S /(j) L 
(Figure^), we find 

a K, s / (k s \ 1 0057 

~%L * ( 1-4413 ± 0 - 0022 ) (*) 

The prefactor 1.4413 implies that smaller grains see more 
of the kinetic stress than larger grains, i.e. they un¬ 
dergo larger velocity fluctuations. Since they are ge¬ 
ometrically less constrained than the larger grains due 




FIG. 3: a: Profile of the kinetic pressure P K (nor¬ 

malised by pgds) for a granular flow with 4> L = 0.30 and 
0 = 23° in steady state; the dotted line shows the best fit 
P K = 0.9Se ~ z ^ 14z ' 8dL ^; b Ratio of the local partial kinetic 
stresses (Jyy S / <Jy]^ as a function of the ratio of the local vol¬ 
ume fractions 4> s /4> L for 4> L = 0.06, 0.15, 0.30, 0.45, 0.60 
and 0.75, and slope 0 = 23°; the plain line shows the best 
fit y= 1.44a: 1 - 0057 

to their size, this result is expected. Interpreting this 
value in the framework of the model by Gray & Thorn¬ 
ton [33], namely introducing a kinetic pressure partition 
coefficient B K such that: 

a y’y L = <t> L 0- + B K <f> s )a^ y , ( 6 ) 

Oyy = (1 - B K (f> L )(T* y , ( 7 ) 

gives a large negative pressure partition coefficient B K ~ 
—0.373, similar to the values observed by m and [39] 
for 3D flows (B k cs —0.39 and B K ~ —0.38 respec¬ 
tively). The value of B K is more then 10 times larger than 
the pressure partition coefficient measured for contact 
stresses in section m and supports the proposition of a 
segregation mechanism driven by granular temperature 
by [SI]. In the chute flow configuration, large grains are 
segregated towards the cooler regions of the flow (namely 
the surface) as predicted by kinetic theory [43] , but differ¬ 
ently from what is observed in rotating drums, for which 
the flow surface coincides with higher granular tempera¬ 
ture [51] . 

We can show that the value of B K is sensitive to the mean 
shear rate. Simulating granular flows with <F L = 0.45, 
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FIG. 4: Kinetic pressure partition coefficient B K as a func¬ 
tion of the flow normalised mean shear rate 7 / yjg/ds com¬ 
puted for flows with a mean volume fraction of large grains 

sag 0.45 and for slope angles varying from 0 — 21 ° 
to 6 — 26°. The dotted line shows the best fit trend 
y = 0.310 + 0.205x. 

and with varying slopes 6 = 21°, 22°, 23°, 24°, 25°, and 
26°, we compute the ratio &yy S /<Jy^ L for each case, and 
fit the dependence on <p s / (j) L by a linear law. We then 
derive the mean value of B K from relations 0 and 0. 
The plot of B k as a function of the normalised shear 
rate is displayed in Figure [4j we observe a monotonous 
increase with 7 / yjg/ds . From the analysis of the segre¬ 
gation time scales for similar 2D numerical flows [32] , this 
increase of B K corresponds to shorter segregation time, 
i.e. a greater segregation velocity. However, it does not 
coincide with a larger segregation rate, as the final posi¬ 
tion of the center of gravity of the large grains was shown 
to remain unaffected by increasing the mean shear rate 

m- 

V. MICRO-STRUCTURAL SIGNATURE 

While one can easily picture why smaller grains are 
submitted to larger kinetic stresses (simply as a result 
of their greater degrees of freedom), the fact that they 
sustain a smaller proportion of the contact stresses than 
their volume proportion is less intuitive. From the micro¬ 
mechanical definition of the contact stress tensor given in 
0. we can try to approximate the partial pressures by 
the mean forces and the mean coordinance number in 
each phase of grains. 

If N C L is the number of contacts involving at least one 
large grain, fs the mean modulus of the forces transmit¬ 
ted by these contacts (in which a least one large grain 
is involved), and is is the mean distance between the 
centres of mass of the two grains in contacts, then an 
estimate U L of the partial pressure P L supported by the 
large grains is (following 0): 

n L ^ 4 v c l /l4. 


FIG. 5: Mean error on the approximations II, II L and II s 
of the pressures P, P L and P s as a function of the volume 
fraction of large grains <j ) L . (For 6 flows with different volume 
fractions of large grains 4> L = 0.06, 0.15, 0.30, 0.45, 0.60 and 
0.75, flowing at 0 — 23°) 

The number of contacts ATf involving at least one large 
grain can be estimated from the number of large grains 
n^, and their coordination number zl (ie the number of 
contacts that a large grain is experiencing on average): 
Nl = n L z L /2. The volume occupied by large grains is 
<p L V] the mean diameter of the large grains being dz,, 
we estimate their number ns = <fi L V/ (ird 2 L / A). Finally, 
£l can be approximated by the diameter of one large 
grain £s ~ d^. This gives the following estimate for 
the magnitude of the partial pressure supported by large 
grains: 

n L ^ \z L h x <7- (8) 

7Td L 

In the same way, being the number of contacts in¬ 
volving at least one small grain, fs the mean modulus of 
the forces transmitted by these contacts (in which a least 
one small grain is involved), and is the mean distance be¬ 
tween the centres of mass of the two grains in contact, 
an estimate of the partial pressure P s supported by 
small grains is: 

n 5 ~ yNsfsis- 

As precedently, Ns = ns x zs/ 2. Smaller grains occu¬ 
pying a volume <p s , their number is ns = (j) S V/(itd 2 s /A). 
Finally, reasoning as for large grains, we suppose £s — ds 
and eventually 

U s ~ \zsh X 0 s . (9) 

7Td S 

Both estimates U s and U L are derived from crude sim¬ 
plifications and are expected to provide only a rough 
guess of the real value of the partial pressures as cal¬ 
culated from 0 . To quantify the error made when using 
^ and ([ 9 ]), we compare the value they provide for P L 
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and P s with the exact computation presented in section 
m averaging over 6 flows with different volume fractions 
of large grains <F L = 0.06, 0.15, 0.30, 0.45, 0.60 and 0.75, 
and flowing at an angle 0 = 23°. The results are shown 
in Figure [5] where | U s — P s | / P s , | II L — P L \ /P L and 
| n-p | jP are reported as a function of </>£. We see that 
in spite of the simplification, II and II L give reasonable 
estimates of the value of P and P L , with a maximum 
error of 3% and 5% respectively, and a mean error of 
2.2% and 1.9% respectively. The error on P s is larger, 
but remains less than 9%, and has a mean value of 6.4%. 
We observe that the error is larger for large proportion 
of mixed contacts between large and small grains, which 
can be seen from Figure [2}-c. As noted before, mixed con¬ 
tacts are the main source of uncertainty when computing 
partial contact stresses. 

If the approximations are too simple to give a quantita¬ 
tive estimate of the stress partition coefficient P, they 
are sufficient to give a qualitative picture of how B be¬ 
haves with the flow micro-structure. Back to the model 
by Gray & Thornton [33], we write 


n L = 4> l (1 + B(j) s ) n, (10) 

u s = c/) s (i — Bcp L ) n, (11) 


as we did for <j£, and <r£, in section 


yy 


yy 


III 


(equations (|3j) 

and ©)• Substituting U L and H s by their expression 
and gives the following for B : 


B ~ 


1 -r d 


zsfs 

ZL f L 


<P L +Td 


Zsfs 

ZlIl 


X (fi‘ 


,S’ 


( 12 ) 


where = d^fds = 2. From this analysis, we see 
that the main ingredient that controls the contact stress 
partition is the ratio of the mean force resultant on small 
grains to the mean force resultant on large grains: the 
smaller this ratio, the more efficient the segregation. For 
rd = 2, the condition for segregation (namely B > 0) is 
that zsf s < \zl / L • We note also that the size ratio 
is not explicitly favorable to stress partition, but tend to 
decrease the value of B. 

This micro-structural condition for segregation can be 
observed from direct measurements of the mean force 
transmitted at contacts, and direct measurement of the 
mean number of contacts in which grains are involved 
on average for each species. Considering a single flow 
(4> L = 0.45, 0 = 23°), we simply count zs and zl the 
number of contacts in which small and large grains are 
involved on average as a function of time (Figure |6]-b). 
For the same flow, the position of the center of mass of 
large grains i/l is reported (Figure [6|a). We see that 
the segregation process implies a significantly larger 
coordinance number for large grains. This is somewhat 
expected since large grains can be surrounded by many 
small neighbours when the converse is not true. The 
difference between zs and zl is maximum when the 
mixing is maximum. More significantly, the mean forces 
fs and /l, transmitted by contacts involving at least 
one small grain and at least one large grain respectively, 
disclose a different behaviour for the two phases of 


a 



FIG. 6 : For a given flow with <f> L = 0.45 and 6 = 23°: a 
Vertical position of the center of mass of the large grains i/l 
normalised by the flow height Ho in the course of time, b 
Coordinance number for zs and zl for small (•) and large 
grains (O) i n the course of time, c Mean force transmitted 
by contact involving at least one small grains (•) and one 
large grain (O) normalised by the weight of one small grain 
msg (fs and /z, respectively, main graph) and normalised 

by the pressure seen by their center of mass (fs and 
respectively, inset) as a function of time. 


grains. The normalised mean forces fs = fs/™s9 and 
f L = fhl^sQ (where ms is the mass of one small grain) 
are plotted in the course of time in Figure [6}c. At the 
start of the flow, the larger grains are at the bottom and 
are thus submitted to a larger pressure than the smaller 
grains closer to the surface: the fact that fs < fL is 
expected. However, as segregation proceeds, and larger 
grains rise to the top, the inequality remains true. This 
effect becomes more apparent when normalising the 
mean forces fs and by the pressure seen by the 
center of mass of each species ( ys and i/l respectively): 

fs = fs/pg(Ho - ys) and f L = f L /pg(H 0 - y L )• Doing 
so, we filter out the influence of the grains position in the 
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flow, and stress the fact that forces transmitted by con¬ 
tacts involving large grains are significantly larger than 
forces transmitted at contacts involving small grains only. 


VI. SUMMARY AND DISCUSSION 

The aim of this work is to understand the size 
segregation mechanism in granular chute flows in terms 
of pressure distribution within the two phases of grains, 
large and small, and to relate it to the granular micro¬ 
structure. Therefore, discrete numerical simulations 
of bi-disperse granular chute flows are systematically 
analysed. Based on the theoretical models by Gray & 
Thornton [33 and Hill & Tan [31 , we compute the 
partial stress tensors associated to the phases of small 
and large grains while separating the contribution of the 
contacts stresses (resulting from the forces transmitted 
at contacts) and the kinetic stresses (resulting from the 
grains velocity fluctuations). Comparing the contact 
pressure in the phases of small grains and large grains, 
we observe a slight deviation from a classical mixture, 
whereby small grains see less of the gravity gradient 
than their volume fraction implies. This result supports 
the possibility of a gravity driven segregation mechanism 
as proposed by [33] . We show however that the contact 
stress partition is extremely sensitive to the definition 
of the partial stress tensors, and more specifically, to 
the way mixed contacts (namely contact involving a 
small grain and a large grain) are handled. As a result, 
the contact stress partition coefficient is sensitive to the 
ratio of grain sizes. This makes the interpretation of 
contact stress ratios awkward, particularly if one wants 
to compare mixtures with different grain sizes: in that 
case, it will be difficult to tell the bias introduced by 
the ponderation by the grains size in the partial stresses 
computation from the influence on the mixed stress 
partition that derives. 

By contrast, the computation of the partial kinetic stress 
tensors is more robust since the separation between the 
phases of large and small grains, relying on the grains 
themselves and not on the contacts, is straightforward. 
Comparing theses partial kinetic stress tensors, we find 


that the phase of small grains exhibits a significantly 
larger kinetic pressure than implied by their volume 
fraction. Computing a kinetic stress partition coefficient 
in the spirit of [33] . we find a value much larger than 
that obtained for contact stresses, and comparable to 
the values found in previous work by [31] and [39!. These 
results support the existence of a segregation mechanism 
induced by shear rate gradients as advocated by [311137], 
and suggest moreover that this mechanism might be 
more important in amplitude than gravity induced 
segregation. In our systems, large grains are segregated 
towards the cooler regions of the flow (the surface) as 
predicted by kinetic theory [43]. In addition, we evidence 
an increase of the value of the kinetic stress partition 
coefficient with the flow shear rate, corresponding to 
smaller segregation time, while no such influence was 
observed on the contact stress partition. 

Finally, using a simple approximation for the contact 
partial stress tensors, we investigate how contact stress 
partition relates with the flow microstructure. We 
observe that grains coordinance (namely the grains 
average number of contacts) is greater in the phase of 
large grains, and that contacts involving at least one 
large grain transmit forces of greater amplitude than the 
others. 

Our results do not allow us to decide for either 
gravity-induced mechanisms or shear gradient induced 
mechanisms, as being the most likely. However, the 
following points should be stressed: i) the computation 
of partial contact stress tensors implies a real difficulty 
due to the existence of mixed contacts, ii) as a result, the 
analysis of partial contact stress ratios lacks robustness, 
iii) the analysis of the microstructure might provide a 
practical proxy to study gravity-driven segregation. 

On the other hand, the robustness of kinetic stress com¬ 
putation, and the sensitivity of the partial kinetic stress 
partition to the flow state, suggest that a parametric 
study of these quantities varying contact properties, 
grains sizes and flow dynamics, would lead to interesting 
new insights into segregation. 
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